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k ship  following  in  the  wako  of  another  ship 
of  comparable  size  experiences  a change  in  hydrodynamic 
force  and  moment  due  to  the  tubulenc©  and  waves  generated 
"by  the  leading  ship.  For  two  identical  wall“3ided  ship3s 
of  finite  draftp  in  tandem?  the  force  and  moment  due  to 
wave 3 f experienced  by  the  trailing  ship*  are  derived 
by  using  bagaily’s  theorem  and  the  usual  assumptions! 
used  5n.  the  theory  of  wave  resistance*  The  force  is 
taken  in  two  components,  horizontal  and  vertical?  and 
the  moment  is  found  about  a transverse  axis  through  the 
of  mass  of  the  trailing  ship* 

Procedures  applicable  to  the  evaluation  of 
all  integrals  involved  are  indicated  and  the  results 
for  the  force  components  and  significant  moment  coxa^ 
ponents  are  listed,  Kijmeria®!  values  are  calculated 
for  som « significant  forcM*  and  moment  components  &nci 
curves  are  plotted  showing  the  trend  of  the  various 
components. 


I*  Introduction 

In  1950*  tha  author a * presented  a paper * 

"The  Dynamic  Stability  on  Course  of  a Towed  Ship"* 
before  the  Society  of  Naval  Architects  and  Marine  Engineers 
(1).  It  was  an  Introductory  Investigation  on  the  stability 
of  a towed  ship  In  which  only  the  motion  in  the  horlsontal 
plane  was  considered . Due  to  the  complexity  of  the 
problem  and  to  determine  the  effect  of  quantities  of 
major  Importance*  many  refinements  were  idealised  or  neg- 
lected* 

One  effeot  which  was  neglected  was  the  hydro- 
dynamic  Influence  of  the  towboat  on  the  towed  ship*  As 
was  pointed  out  by  £•  H,  Peters*  In  the  discussion  Which 
followed  the  afore  mentioned  paper*  this  effeot  Is  very 
significant  In  the  analysis  of  ths  towing  problem, 
especially  for  the  case  of  a short  tow.  This  fact  was 
experimentally  verified  by  B.  0.  Barrlllon  (2)  whe 
found  that  the  hydrodynamic  Influence  on  the  trailing 
ship  was  noticeable  even  when  the  distance  between  the 
towing  and  towed  ship  was  of  the  order  of  five  ship 
lengths. 

To  obtain  a physical  pioture  of  ths  problem* 
imagine  an  unpowered  ship  moving  with  steady  velocity 
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on  s.  calm  sea*  The  ship  experiences  & total  resistance, 
the  sum  of  frictional  find  residual  resistance,  of  which 
«a?s  resistance  is  the  chief  part  of  the  latter*  If 
now  another  unpowered  ship  mowing  with  the  same  steady 
velocity  is  introduced  directly  ahead  of  the  existing 
ship,  the  trailing  ship  will  experience  a change  in 
resistance*  The  Change  in  resistance  is  largely  due 
to  the  turbulenee  end  waves  generated  by  the  leading 
ship*  The  turbulence  and  additional  waves  affect  both 
the  frictional  and  wave  resistance  encountered  by  the 
trailing  ship.  In  addition,  the  presence  of  the  leading 
ship  near  the  trailing  ship  will  give  rise  to  forces  of 
mutual  action  which  may  be  classed  as  dive  to  local 
disturbances • 

In  1936,  T,  H.  Havelock  wrote  a paper  on  the 
mutual  aotion  between  two  bodies  (3)*  The  general  aotion 
between  two  bodies  are  given  and  are  applied  to  several 
oases  of  a pair  of  doublets  which  are  oriented  in 
different  positions*  Only  the  horisontal  component  force 
; is  do alt  with  and  the  vertical  component  and  moment 
experienced  by  th»  trailing  body  are  not  discussed^ 

In  a later  paper  (2),  Havelock  applied  one 
part  of  the  equations  developed  in  hie  earlier  paper  to 
a wall-sided  vessel  of  infinite  draft*  Equations  for 
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the  horizontal  component  forces  wore  given  but  curves 
showing  the  trend  were  not  sho^m.  As  a first  approxi- 
mation, the  use  of  an  infinite  draft  is  quite  adequate 
since  those  parts  of  the  shJp-form  which  are  situated 
very  deep  below  the  fluid  surface  do  not  affect  the 
wave  resistance  seriously*  For  a closer  approximation, 
however,  Havelock  does  show  (4)  that  the  effect  of 
decreasing  draft  becomes  appreciable  especially  for 
email  draft-length  ratios  at  high  velocities* 

The  present  investigation  is  an  extension  of 
the  work  done  by  Havelock  and  deals  with  only  that  part 
of  the  hydrodynamic  force  and  moment  experienced  by  the 
trailing  ship  which  is  due  to  waves  and  the  presence 
of  the  leading  ship*  The  finite  draft  which  was  neglected 
by  Havelock  will  also  be  considered*  The  effect  of 
fr5.ctlon  on  the  force  and  moment,  which  Is  complicated 
by  the  turbulence  created  by  the  leading  ship,  will  not 
be  discussed* 
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II.  Assumptions  and  Methods  in  Wave  Resistance 

Any  Burfaoe  ship  In  motion  encounters  a total 
resistance  whioh  m&„  be  classified  Into  several  mutually 
dependent  components.  They  may  be  classified  as  (a) 
frictional  resistance  which  is  assumed  to  be  the  resistance 
of  a flat  plate,  of  the  same  roughness,  length,  and 
surface  area  as  the  hull  of  the  ship;  (b)  pressure 
resistance  caused  by  pressure  changes  due  to  viscosity, 
especially  in  the  after  body;  and  (c)  wave  resistance 
which  follow  behind  a ship.  Another  type  of  resistance 
is  spray  resistance.  However,  this  is  an  independent 
physioal  effeot  and  is  only  Important  for  high  speed 
ships* 

In  the  present  state  of  development  the  the* 
oretlcal  analysis  of  wave  resistance  is,  at  best,  an 
approximation.  To  facilitate  practical  analysis  many 
assumptions  have  to  be  made  and  thus  a complete  solution 
cannot  be  obtained  even  for  the  simplest  of  ship-shaped 
forms.  Perhaps  the  greatest  violation  of  physical  fact 
is  the  use  of  an  ideal  fluid  end  the  assumption  that 
the  wave  resistance  is  independent  of  the  other  forms 

» 

of  resistance.  These  assumptions  are  justifiable  only 
because  of  the  need  for  systematically  simplifying  a 
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complex  problem,. 

Further  assumptions  aro  as  follows: 

(1)  The  velocity  deviations  due  to  the  presence 
of  the  ship’s  hull  are  small  relative  to 
the  main  stream  velocity, 

(2)  The  height  of  waves  generated  by  the  3hip’s 
hull  is  small  relative  to  their  length, 

(3)  Trim  and  linkage  do  not  affect  the  wave 
resistance. 

(4)  Interference  between  the  ship’s  hull  and. 
wave  system  does  not  affect  the  wave 
resistance. 

(5)  Tbs  inclination  of  the  tangent  plane  to 
tho  median  plane  et  any  point  on  the 
surface  of  the  ship  is  small, 

(6)  The  fluid  Is  non«viaeoue„  of  infinite 
depth  and  unbounded  on  the  free  surface,, 

(7)  Surface  tension  on  the  free  surface  is 
negligible, 

(8)  The  condition  of  atmospheric  pressure  on 
the  free  surface  Is  satisfied  at  the 
undisturbed  free  surface  and  not  at  the 
disturbed  free  surface. 


«*»#t** 
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Takiag  all  these  assumptions  Into  consideration* 
the  actual  physical  picture  is  far  removed  from  the 
theoretical  picture  with  its  many  assumptions*  With 
all  these  assumptions*  the  general  agreement  which  exists 
between  experimental  and  theoretical  results  is  very 
surprising  indeed* 

There  are  several  methods  of  calculating  the 
effect  of  waves  on  the  hull  of  a shipe  To  mention  a 
few*  the  pressure  changes  due  to  wave  motion  may  be 
integrated  over  the  surface  of  the  hull}  the  wave  energy 
left  behind  the  ship  may  be  calculated;  the  energy 
dissipation  using  Rayleigh's  dissipation  coefficient  may 
be  calculated;  and  forces  experienced  by  hydrodynamic 
singularities  (sources  and  sinks)  may  be  calculated*  , 

The  method  which  seems  best  suited  to  the  two  body 
problem  of  towing  Is  the  last  one*  the  calculation  of 
forces  experienced  by  the  hydrodynamic  singularities 
which  together  with  the  free  stream  represent  the  bodies* 

Legally* s theorem  (5)  is  indispensable  in  this 
method*  It  enables  the  calculation  of  forces  and  moments 
experienced  by  the  hydrodynamic  singularities*  The 
resultant  forces  Pi  and  resultant  moments  Mi  acting 
on  a body  whose  surface  is  a closed  stream  surface  of 
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the  fluid  mpt:lon,  are  given  by  the  vectors 


% = - 4TTp  miqi 
Mi  a -t*i  X 4 T T p niqi 


(2.1) 


where  mi  is  the  strength  of  a source  internal  to  the 
stream  surfaoe,  q^  is  the  resultant  fluid  velocity 
vector  at  the  location  of  the  source  due  to  all  other 
souroes,  p is  the  mass  density  of  the  fluid,  and 
Fi  is  the  displacement  of  the  souroe  aii  from  the 
moment  center.  The  forces  Fi  have  the  direction  of 
- qi  and  the  lines  of  action  pass  through  the  point  at 
which  the  source  is  located.  The  vectors  Mi  which 
represent  the  resultant  moments  is  perpendicular  to  the 
plane  formed  by  r|  and 

The  motion  of  the  ideal  fluid  produced  by 
the  hydrodynamic  singul arities  and  the  free  stream^s 
described  by  a velocity  potential  ^ (x,  y,  z,  t)  which 
must  satisfy  Laplace's  equation  and  Conditions 

at  the  boundaries  of  the  fluid,  at  the  surface  of  the 
hull,  and  at  the  free  surface  of  the  liquid.  It  must 
also  satisfy  the  aquations  of  wave  motion.  Of  all  the 
conditions  to  be  satisfied,  the  conditions  at  the  free 
surface  cause  the  most  difficulty. 
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III,  Free  Surface  Condition 

In  elementary  hydrodynamics  problems,  any 
body  moving  uniformly  in  an  unbounded  ideal  finis  at 
rest  at  infinity,  experiences,  by  D* Alembert's  paradox, 
no  resistance.  It  may,  due  to  its  orientation  with  the 
free  stream,  experience  a moment.  Unlike  an  unbounded 
fluid,  in  the  practical  application  of  theoretical 
hydrodynamics  to  problems  of  naval  architecture,  a free 
surface  must  always  be  considered.  For  this  case, 
despite  the  use  of  an  Idealised  fluid  in  the  analysis, 
a body  iuovlng  uniformly  on  or  near  a free  surface  en- 
counters resistance  due  to  the  energy  required  in  creating 
gravity  waves  which  travel  with  tho  body. 

To  obtain  the  free  surface  condition,  consider 
a stationary  right  handed  systen  of  xys  axes  in  the  free 
surface  of  a stream  flowing  with  steady  velocity  o in 
the  negative  x-direction.  The  origin  and  the  x and  y 
axes  are  plaoed  in  the  undisturbed  free  surface  and  the 
z-axis  ia  vertically  upward. 

Now  consider  two  points  on  the  free  surface  of 
the  stream*  One  point  is  taken  near  a disturbance  and 
the  other  point  is  taken  at  Infinity.  The  disturbance 
ia  of  such  a nature  that  its  effect  at  infinity  is 
negligible.  For  these  two  points,  assuming  that  any 


deviation  of  a fluid  particle  Is  resisted  by  a force 
proportional  to  the  fluid  velocity,  the  pressure  equation 
may  be  written  as 

7"  +f-  (3.1) 

where  p Is  the  mass  density  of  the  fluid,  p ana 

p0  are  the  pressureo  at  the  disturbance  and  at  Infinity 

respectively  and  >{  is  the  surface  elevation  from  the 

undisturbed  free  surface.  The  quantities  - (c  -t*  <f>  ), 

- 4>u,  are  the  absolute  velocities  of  the  fluid  at 

•j  12 

the  disturbance.  The  quantity  ft' $ is  the  force  poten- 
tial which  introduces  the  hypothesis  of  a frictional 
force  proportional  to  the  deviation  of  the  fluid  velocity 
from  the  uniform  flow  c.  This  hypothesis  is  adopted  to 
keep  various  integrals  convergent  in  the  subsequent  analy- 
sis (6),  (7),  However,  since  wave  resistance  is  calcu- 
lated for  a friotionless  fluid  the  final  result  will  be 
simplified  by  letting  the  frictional  coefficient  /t'  tend 
to  zero. 

The  quantities  -{»  , “4U»  “4-  a**®  the  devia- 

tions,  at  the  disturbance,  of  the  fluid  velocity  from 
the  uniform  flow  c in  the  x,  y and  z directions. 

By  assuming  that  the  vector  sum  of  these  deviations  Is 
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small,  l.o,,  that  ( + $>*)  is  negligible,  the 

pressure  equation  becomes, 

T+C^  + f)?V4’=^-  (3.8) 

Now  by  assuming  that  the  pressure  on  the  free  surface  at 
the  disturbance  is  equal  to  the  pressure  at  the  undis- 
turbed surface  at  infinity,  i.a.,  p » pQ  , 


c4>„  + g >1  “/*'<£  = O 
or 

^ c^x)  i3*3> 

This  equation  gives  the  elevation  of  the  free  surface  at 
the  disturbance. 

The  general  boundary  condition  wi  ich  must  be 
satisfied  at  the  free  surface  is  that  there  shall  be  no 
flow  across  the  free  surface,  i.e.,  the  velooity,  rela- 
tive to  the  free  surface,  of  all  particles  lying  in  it 

U M !•  A>  a — , i»_  - - V 4k  i.u  - A # _ 

muou  w«U56uw  vu  uao  ouTi  11  UUQ  1 I OU  SUTJL  HUO  J.S 

given  by  the  equation  F(x„y,z,t)  = 0,  this  boundary 
condition  may  be  expressed  as 


— = F.  +•  uF,  vF_  + wF*  » 

at  * x y 
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where  u - (c  -r^),  v w re  “ 4^* 

In  the  present  problem,  the  equation  of  the 
froe  surface  is  given  by 

* « ^ (x,  j) 

or 

F(x,y,z,t)  « z — >j  U,y).  (3.5) 

The  Eq.  (3.4)  becomes 

<3-6) 

Assuming  that  the  wave  height  is  small  in  comparison  to 
the  wave  length,  the  slope  of  the  waves  become  small 
quantities  of  first.  Then  by  neglecting  small  quantities 
of  second  order,  Eq.  (3.6)  becomes, 

ix‘  <3-7> 

Combining  Eq.  (3.3)  and  Eq.  (3.7)  the  result  is  the  free 
surface  condition 

+w+  *»\ -/*♦„  *0  (3-8) 


r~ 


at  z =»  0.  and  whore  g/c^  and  =At'/c.  The 

velocity  potential  which  describes  the  motion  of  the 
ideal  fluid  and  the  distribution  of  sources  and  sinks, 
which  together  with  the  free  stream  represents  the  bodies, 
must  satisfy  this  free  surface  condition. 
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IV.  Velocity  Potential 

Using  the  same  axes  as  in  the  previous  section, 
consider  a siaple  source  of  strength  m at  (0,0, -f), 
traveling  with  uniform  velocity  c in  the  positive 
’■-direction.  In  the  absence  of  a free  surface  and  in 
a fluid  of  infinite  extent  in  all  directions,  the  velo- 
city potential  is 

<|>  - (4.1) 

where  (r^)*5  * x2  y2  + (*  -t  f)2.  In  the  presence  of 
a free  surface,  the  velocity  potential  may  be  writ ton 
as  (8) 

4>  » ~'/rl  + ^ “i/1*!  (4.2) 

where  (rA)2«s  (x  - x^)2  *+■  (y  - y^)2  -♦-(*-  s^.)2. 

The  summation  ^ m«/r<  represents  an  image  system, 

\ 1 1 

a distribution  of  hydrodynamlo  singularities,  in  ar> 
unbounded  fluid,  which  together  with  the  given  souroe 
will  produce  a free  surface. 

Making  uuie  of  the  known  integral  (8) 
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(x2  + y2  + z2f*  ™ 
for  z > 0 and 
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the  velocity  potential  is  assumed  to  be  of  the  form 
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where  55  * x cos  0 t y sin  0,  and  only  the  real  part  is 
to  be  used. 

Substituting  <p  into  the  free  surface  condi- 
tion, Eq.  (3.3),  it  can  be  shown  that 


F<e,K)  = 


rp 

2TT 


r z / K +iM*cce  \ 1 
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(4.5) 
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The  velocity  potential  than  becomes,,  after  soma  reduc- 
tion. 


*2  _ £»  K,  r*i 
‘r,  fi  ” TT 


(•"fc-fHisJ 

j " 

K-»f0secld  + i>iAec© 
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f 4.6) 


where 


to  •= 


x2  -r  3T2  Hr  (a  -j-  f)2 
x2  + y2  + (z  - f)2 
x cob  0 + y sin  0, 


The  limiting  value  of  Eq»  (4-6)  la  to  be  taken  for 

The  first  term  represents  the  given  source  at 
(OsO-f);  the  second  term  represents  a sink  at  (0„0,f), 
and  the  last  term  represents  a continuous  distribution  of 
sources  and  sinks,  in  the  plane  z =*  f , trailing  in  the 
negative  x-direction  to  infinity,  The  last  two  terms 
form  the  image  system  which  is  required  to  produce  the 
froe  surface  (9). 


Instead  of  an  isolated  5uuiC6,  f or  a general 
continuous  distribution  of  sources  below  the  free  sur- 
face, the  source  strength  m is  replaced  by  J*  <r  dS 
.-r—  where  «j*  is  the  surface  density  of  a source  on  a surface 
dS  at  a point  (h,k,-f)  within  the  fluid.  Applying  the 
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principle  of  Inverse  flow,  the  velocity  potential  for  a 
general  continuous  distribution  of  sources  in  an  ideal 
fluid  with  a free  surface  is 
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whore  o is  the  speed  of  the  free  stream  in  the  negative 
x-direction  and 


rj,2*:  (x  - h)2  -r  (y  - k)2  + (x  -f-  f)2 

r22  «.  (x  - h)2  r (y  - k)2  -r  U - f )2 

£5  o (x  - h)  cos  0 -t*  (y  - k)  sin  0. 


«*&*■  'imrr7.r*fF  raw  1 > 


( V.  Source  Distribution 

A complete  determination  of  the  distribution 
of  souroes  and  sinks  which  will  satisfy  all  the  boundary 
conditions  is  a very  difficult  problem*  The  problem 
is  somewhat  simplified  by  considering  a distribution  of 
sources  and  sinks  taken  over  the  plane  y « 0 instead  of 
a distribution  over  a volume*  This  is  only  permissible 
for  a narrow  ship’s  hull*  If  at  the  point  (h*0*-f)»  <r 
is  the  source  density*  the  total  atx’ouguu  ov*?  a csidl 
area  dS  is  (TdS  and  the  corresponding  flux  is  4TTC“dS, 
Considering  the  entire  souroe  distribution*  since  it  is 
continuous  within  its  limits*  the  velocity  flows  outward 
normally  on  both  sides  of  the  plane  y « O-  The  velocity 
normal  to  the  plane  ys  0 i»  thus  4 IT  0“  dS/  2dS  « ± 2 T (T  • 
This  normal  velocity  is  related  to  the  form 
of  the  ship  by  considering  the  condition  which  must  be 
satisfied  at  the  surfaoe  of  the  ship’s  hull*  The 
boundary  condition  is  that  the  normal  velocity  of  the 
ship’s  hull  and  the  fluid  in  oontace  with  it  must  be  the 
same*  The  equation  which  expresses  this  boundary  condition 
is  the  same  as  Eq*  (3*4)* 

pt  + ♦ TPy  +*  wFx  * 0 
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If  y » F(x,z)  Is  the  equation  of  the  surface  of; the 
ship,  then  this  boundary  condition  is  expressed  by 

(o  - u)  Px  + v • w P8  * 0 (5,1) 

It  is  assumed  here  that  the  inclination  of  the 
tangent  plane  to  the  median  plane  at  any  point  on  the 
surface  of  the  ship  is  small.  This  condition  is  not 
satisfied  at  the  bottom  of  modern  ships.  However,  since 
most  of  the  wave  resistance  is  due  to  the  upper  part  of 
the  immersed  portion  of  the  hull,  it  is  a reasonable 
assumption.  With  this  assumption  Fx  and  Fz  will 
be  small,  and  since  u « - and  w ~ ~ 4* 

were  considered  small;  by  neglecting  small  quantities  of 
second  order,  Eq.  (5.1)  becomes 

v * - e Fx.  (5<>2) 

How  by  assuming  that  the  velocity  in  the  y»direction  at 
the  hull  is  equal  to  the  velocity  normal  to  the  plane  y —0, 


(T 


Z7T  * zTT  ^X 


(5.3) 


where  dy/dx  is  obtained  from  y F(x,s)  the  equation 
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describing  the  surface  of  the  hull.  Eq.  (5.3)  Is  then 
the  usual  approximation  for  the  density  of  the  distribu- 
tion of  sources  and  sinks  over  the  median  plane  of  the 
ship’s  hull. 
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VI o Forces  cn  the  Towed  Ship 

Consider  two  identical  wall-sided  ships  A 
and  B,  oi. r right  behind  the  other,  placed  in  a stream 
which  flows  with  constant  speed  c in  the  negative 
x-direction.  The  ships  are  spaced  such  that  the  distance 
between  their  centers  of  masses  is  equal  to  L,  l'he  ships 
are  of  length  2 l , beam  2b,  and  draft  d.  Ship  A rep- 
resents the  towed  ship  and  ship  B represents  the  towing 
ship. 

A right  handed  system  of  body  axes  x,.y*  z 
with  the  origin  In  the  undisturbed  free  surface  is 
fixed  in  ship  A as  shown  in  Fig.  (5.1).  The  equation 
which  describes  the  water  plane  section  of  ship  A is 

7 - b [l  - (x/l  )2]  (6.1) 

where  i x * l . For  the  ship  B the  equation  is 

yeb^l-  £(x-l»)/lj2|  (6.2) 

where  (I*  - l ) ^ x ^ (L  + l ),. 

Each  ship  will  now  be  replaced  b a continuous 
distribution  of  sources  and  sinks  over  the  vertical  median 
plane  as  described  in  the  last  section.  Following 
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Eq.  (5*3)  and  referring  all  quantities  to  the  axes  fixed 
in  ship  A,  the  source  densities,  crA  and  c , for  the 
ships  A and  B are 

<71  & (ho x)/(7T£2) 

A 

where  -l  £ x £ l , and 

= [bc(x  *•  L)j  /(7TZ2) 

where  (L-Z)^X^(L-t-Z)* 

Prom  Eq,  (4,7)  and  again  referring  all  quantities 
to  the  axes  fixed  in  ship  A,  the  veloolty  potential  of  the 
above  source  distribution  in  the  uniform  stream  will  be. 
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where  u and  w are  the  magnitudes'  of  the  components  of 
<f  In  the  x and  % directions  respectively*  The  result 
tant  forces  In  the  x and  & directions  on  the  distrihu~ 
tlon  A will  then  be 
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o“gdh*  df  * • 
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The  velocity  components  u and  w in  terms 
of  the  velocity  potential  $ are 


(6.8) 


where  the  term  (1/r^)  Q'^dhdf,  the  velooity  potential  for 
a source  at  any  point  fx,y9s),  has  been  eliminated  from 
• Then 
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5 B f ttJfi - , 

[(h»-  h)2  + (r  -»*  f»)e]% 

J = X ■..(! . 

Z [(h»-  h)2  + (f  - f*)2]V* 


h (6«u) 


Substituting  Eq,  (6.9)  and  Eq,  (6*10)  into  Eq,  (6*7) 


Rx  s X1  + ^ XS  * X4  X5  * V 


(6,12) 


where 
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The  result  for  the  force  R8  Is 


Rz  = Z1  ^ zg  +■  z3  + z4  ■*"  z5  * 


(6*13) 
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VII.  Interpretation  of  Force  Components 

Tho  different  components  of  the  forces  In  the 
x and  z directions  represent  the  forces  between  the 
various  parts  of  the  source  distribution.  To  describe 
the  various  force  and  moment  components , the  Interpreter 
tlon  of  the  various  parts  of  the  velocity  potential  which 
describes  the  flow  due  to  the  source  distribution,  which 
is  given  in  Section  IV.  must  be  known. 

?h8  first  component  of  Rx  , 


TTp 


o'* 


(T/Jh'jf' 


represents  the  force  on  the  distribution  A due  to  the 
uniform  stream.  However,  since  the  total  source  strength 
over  the  distribution  is  zero,  it  is  equal  to  zero. 

For  the  second  component 
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consider  two  sources  of  strength  m sncl  mr  (it  the 
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points  (h,0,f ) and  (h*»0,=f*)  respectively*  as 
shown  in  Fig*  (7,1).  The  magnitude  of  the  force  between 
m and  m*  Is  given  by 

4t7> m m* 
r® 

where  r la  the  distance  between  m and  a* . The  force 
acts  along  the  line  drawn  between  m and  m*  and  is  an 
attraction  when  m and  m*  are  of  like  sign.  When  m 
and  m*  are  of  unlike  sign  the  for  >e  Is  a repulsion. 
Referring  to  Fig,  (7.1) 

r2  ~ (h*~  h)2  -r  (f  + f1)2. 

Then  the  force  between  m and  nt'  becomes 

4 TT  f>  m m* 

(h»-  h)2  +■  (f  + f • )2 

The  horizontal  component  of  this  force  is 

AIT  pm  m 1 (h*»  h) 

[(h*-  h)2  +•  (f  + f*)2]^ 


I i*  a 


represents  a certain  distribution  over 


an  area  A in  the  plane  y « 0,  we  may  writs 


<rxdhdf 


Similarly  if  a*  represents  an  identloal  distribution 
except  possibly  for  the  sign  of  ra* y we  may  write 


m* 


<TA*dh*df*. 


Then  for  the  given  distribution  of  sources,  the  magni- 
tude of  the  resultant  force  in  the  x-dlrection  is 


(TjJJ!  if 


Comparing  this  equation  with  the  component  it  is 

clear  that  the  form  and  magnitude  are  the  same.  Hence 
Xg  represents  the  resultant  force  in  the  x-dlrection, 
on  the  distribution  A,  due  to  the  forces  which  exist  be- 
twoen  the  sources  and  sinks  in  the  distribution  A and 
its  image  system,  excluding  that  part  of  the  image  sys- 
tem which  trails  aft  to  infinity. 
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By  following  the  seme  reasoning.  It  can  be 
shown  that  the  component  X3  la  the  resultant  force 
In  the  x-dlrectlon,  on  the  distribution  A,  due  to  the 
forces  between  the  sources  and  sinks  In  the  distribution 
A and  the  Image  of  the  sources  and  sinks  In  the  distri- 
bution B,  again  excluding  that  part  of  the  image  system 
of  the  distribution  B which  trails  aft  to  infinity.  The 
component  X^  represents  the  resultant  force  in  the  x- 
dlrectlon,  on  the  distribution  A,  due  to  the  forces  be- 
tween the  sources  and  sinks  in  the  distribution  A and  B» 
For  the  remaining  distribution  which  trails  aft 
to  infinity  for  both  distribution  A and  B,  the  velocity 
potential  is 


r pdpl  1 pT  P*-,«f+miKS 

— <n<lhJf+  CTlofJf  S€c"6 d&  -£ — 

T A a 9d  + i/AS< 

J 4.1  *^o 


as  given  by  Eq*  (6.4).  The  first  part  gives  the  trailing 
system  for  the  distribution  A and  the  second  part  gives 
the  trailing  system  for  the  distribution  B. 

The  resultant  force  in  the  x-dlrectlcn  on  the 
distribution  A for  this  case  cannot  be  shown  as  simply 
as  was  done  in  the  previoua  oases.  However,  a little 
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reflection  on  Legally’ s theorem  and  this  part  of  she  velo- 
city potential  will  show  that  Xg  represents  the  resul- 
tant foroe  in  the  x-direetion,  on  the  distribution  A , due 
to  the  forces  between  the  trailing  system  for  the  distri- 
bution A and  the  distribution  A itself.  This  expression 
is  the  usual  form  of  the  wave  resistance  for  a ship  alone 
on  a calm  sea. 

The  component  Xg  is  the  resultant  force  in 
the  x-direction  due  to  the  forces  between  the  distribu- 
tion A and  the  trailing  system  for  the  distribution  B. 
Physically  this  term  represents  the  wave  resistance  on 
the  distribution  A duo  to  the  wave  interference  caused 
by  the  waves  left  by  the  distribution  B,  or  the  towing 
ship. 

With  the  foregoing  interpretation,  the  force 
Rx  acting  on  the  distribution  A,  may  be  classed  into 
three  parts;  (1)  the  wave  resistance  of  the  towed  ship 
a8  if  existing  alone  on  s calm  sea,  (2)  the  mutual 
action  between  the  two  systems  which  may  be  classed  as 
due  to  local  disturbances,  and  (3)  the  wave  interference 
acting  on  the  towed  ship  (3). 

Following  the  same  procedure,  the  interpretation 
of  the  various  components  in  Ra  may  be  given.  The  first 
component  represents  the  resultant  vertical  force  on 


the  distribution  A due  to  the  forces  between  the  distri- 
bution A end  its  image  system  excluding  that  part  of  the 
image  system  which  trails  aft  to  infinity.  The  compo- 
nent Zg  represents  the  resultant  vertical  force  on  the 
distribution  A due  to  the  forces  between  the  distribu- 
tion A and  the  image  system  of  the  distribution  B,  again 
excluding  that  part  of  the  image  system  of  the  distribu- 
tion B which  trails  aft  to  infinity.  The  component  Z3 
represents  the  resultant  vertical  force  due  to  the  forces 
between  the  distribution  A and  the  distribution  B.  Finally, 
Z4  and  Z5  represent  the  resultant  vertical  force  on  the 
distribution  A due  to  the  forces  between  the  trailing 
system  for  both  distribution  A and  B and  the  distribution 
A itself. 
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VIII*  Moment  on  Towed  Ship 

To  obtain  the  expression  for  the  moment* 
consider  again  the  source  of  strength  o-A*dh’ df1  at 
the  point  (h*  *0*~f»)  in  the  distribution  A,  By 
Legally7 a theorem  for  momenta*  the  moment  about  the 
center  of  mass  of  the  distribution  A due  to  the  force 
on  this  source  la  given  by 

dL«-rX4T»p(  crA»dh*df* ) q (8.1) 

where  r is  the  displacement  veotor  of  the  source 
CT^'dh'df*  from  the  center  of  mass  of  the  towed  ship  as 
shown  in  Fig.  (8.1). 

If  the  source  distribution  is  taken  in  the 
vertical  median  plane  of  the  towed  ship 


r » h*  i + (f  * - e ) k 


(8.2) 


"i  and  k are  unit  vectors  in  the  positive  x 
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of  the  center  of  gravity  of  the  towod  ship  from  the 
origin  of  the  coordinate  system.  The  resultant  velocity 
at  the  point  (h^O.-'f*}  may  be  written  as 


q = u 1 t y J t t k 


(8.3) 
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MASS  CENTER  AND  DISPLACEMENT 
VECTOR  r FOR  MOMENT  . 
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where  j is  a unit  vector  in 
for  a distribution  of  sources 
plane i,  the  resultant  velocity 
equal  to  sero.  Then 


the  y-direction.  However* 
taken  in  the  vertical  median 
v at  this  point  will  be 


r X q «■ 
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and 
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(8.5) 


The  total  resultant  moment  on  the  distribution  A about  an 
axis  through  the  center  of  gravity  and  parallel  to  the  y- 
axis  will  then  be 
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The  quantities  u and  w are  given  by  Eq.  (6.9) 
and  Eq.  (6,10)  respectively.  Substituting  into  Eq.  (8.6) 
the  moment  then  becomes 


Ly  - ••••  + 
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where 
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As  was  dons  for  the  forces  Rx  and  Ra  the 
interpretation  for  the  various  components  of  is  as 
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follows.  represents  the  moment  on  the  distribution  A 

about  the  center  of  mass  due  to  the  uniform  stream*  The 
sum  of  the  components  Mg  through  Mg  represents  the 
moment  of  the  distribution  A about  the  center  of  mass  due 
to  the  horizontal  components  of  the  forces  between  the 
distribution  A and  (a)  the  image  system  of  the  dlatribu- 
tion  A excluding  the  trailing  system  of  sources  and  sinks, 
(b)  the  distribution  B»  (o ) the  image  system  of  sources 
and  sinks,  (d)  the  trailing  system  of  the  distribution  A, 
and  (e)  the  trailing  system  of  the  distribution  B.  In 
like  manner  the  sum  of  the  components  M7  through  Mn 
represents  the  moment  due  to  the  vertical  component  of 
the  forces  between  the  various  parts  of  the  source  dis- 
tribution. 
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IX*  Reduction  of  Some  Integrals 

Many  of  the  Integrals  which  appear  in  the  ex- 
pressions for  Rz  and  may  be  rewritten  In 

.other  forms  more  suitable  for  further  Integration*  In 
all  of  these  expressions  we  have  the  Integrals 
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These  integrals  may  be  evaluated  by  considering 
K as  a complex  variable  and  using  a suitable  contour* 
Both  5ntegrals  will  have  simple  poles  within  the  oontours 
and  the  position  of  the  simple  poles  will  depend  upon  the 
sign  of  see  9.. 
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pole  of  the  integral , the  range  of  integration  with 
respect  to  0 Is  reduced  from  -IT  , T to  0,  fT /2  » 
Then  the  Integration  In  0 and  K may  be  written  as 
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Choosing  a contour  bounded  by  the  positive 
half  of  the  real  axis  and  the  negative  half  of  the  imagi- 
nary axis,  it  can  be  shown  that 


W f,(6^K)e  KdK  = 


- 2U  i {k.s«*0  j-«iK.lfc*-K)«ee  - K.(# 


« — 


r 


■i*  I 


w (fc-b)<osd 

>t(ht)e  mdm 


, mi  vr  9 ■ w «*<*rrv* i 1 
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where  (h'-  h)  < 0.  Here,  it  Is  to  be  noted  that  /l-*0. 
The  expressions  for  | («)  and  ^ (m)  are  as  follows: 

'l  * 


?» 


ICo^ec*#  [cos  m + ms»M  m (f+'F) 


(9.6) 


ftt  *Cfc*6  fc»in  mCfelOj-rn  gpg  Wlfft-f*) 


Hn*  ♦*  1^*  SttC4d 


For  (h»-  h)  > a Contour  bounded  by  the  posi- 
tive half  of  both  the  real  and  Imaginary  axes  la  ohosen. 
Then 


M-+0 


J *, <»<*)« 


KJK 


>,(ni)e  mim 


(9.7) 


r 

l **(")« 


tnim 


where  again,  /U-yQ* 


rmr^O£s 
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Poll  owing  similar  procedures*  it  may  be  shown 


that 


Lk  tf  ( *,K) 

M-+o  *■ 

J A 


ZTTi  Jc 


-k/MV 

e KdK 


ft,  V m^-»0c®s6 

- ^(w/e  **»4i 

•*© 


-h)  cos  6 


where  (h*-  h)  <.  0,  and 


r«° 

«|»  ($*K) 
.U-*o  1 


= iv> 


~rr»(H,-l» ) f.OS  ® 


. r % -m(Mccs9 

*t  l £j,(m)  e 


II  /* 


where  (hf-  h)  > 0, 
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Substituting  these  equations  into  Eq.  (9.3) 
tn*  result  is 


N1  " 2 [ Sic 9 40  ^a<^co»[^(l-0««ce] 


. P . . rr,(h-h)casf>  , 7 

41i  £»(m)e 


l(9o!0) 


for  (h*-  h)  < 0,  and 


(9.11) 


for  (h‘-  h)  > 0. 

For  Ng  the  result  Is 


N 


2 


} | S4 t*G  J B ^4  TT  c 0 e ^ « jj^  t] 


riv olz) 


for  (h* 


h)  < 0,  find 
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c 

V 


*8= 


sec  6d0 


uJ 


» * / \ "m  th-H)  cos  6 . 

2 £,(*»?«  mditi 


(9*13) 


for  (h»-  h)  > 0 . 

Returning  to  the  expressions  for  Xg  and 
Xg  , It  is  to  be  noted  that  they  are  imaginary.  By 
substituting  the  imaginary  part  of  Eq.  (9.10)  and 
Eq.  (9.11}  into  Xg  and  Xg  , real  forces  are  obtained. 
The  result  is. 


H + X6 


€ 


■ - i L p kJ  Ic^Jh'df'l  I^JdWff  »etOdfl  *'^C*  ^ 


I J 

'o'!  ** 


mdm 


+/fcp*. 


M 


. . . 

5«&Sd6  i»Jm 


tKXplf  sec*  itCC ioafctk'+lsectf  id 

•o'**  **s  » w 


(9.14) 


+ lfcf>K.  I \rM\v>«  *ec6l&  f 5J,(rr»)eW’^4,>eA,emir»» 
6*-J  J* 

-52  Xp  ii'if j**4  * ® e fccw |i^( k'H») ttcfjj J9. 
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Further  reduction  on  this  equation  is  possible. 
Consider  the  first  two  terms  in  the  right  member.  In 
the  first  term,  the  integration  with  respect  to  the  vari- 
ables h and  h’  is  over  the  area  shown  in  Fig.  (9.1). 
By  changing  the  order  of  Integration,  the  Integration  in 
h and  h'  becomes 


>2  p* 
•it  \ 


cot  6 

a 


(9.15) 


By  interchanging  h and  hf,  the  value  of  the  Integral 
la  not  changed.  Hence  the  integration  in  h and  h’ 
may  be  written  as 


/ 


<rAlh 

ti 


m ( h'-b)  e 


(9.16) 


In  expression  (9.15),  there  is  a limitation  on  (h;-  h), 
>i.e.,  (h’-  h)  > 0.  By  the  process  of  interchanging  h 
and  h* , the  limitation  on  (h*«  h)  in  expression  (9.16) 
becomes  (h*-  h)  < 0. 

Now  by  comparing  the  integration  in  h and 
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C 


c 


u 


hT  in  the  second  term  of  the  right  member  of  Eq.  (9.14) 
with  expression  (9.16),  it  is  to  be  noted  that  they  ere 

td.nf.tnal  avnanf  fm>  f.hatn  mitm.  Itina  fcha  aim  nf  fVm 

first  two  terms  in  the  right  member  of  Eq.  (9,14)  is 
equal  to  zero. 

By  s similar  procedure,  the  third  term  in  the 
right  member  of  Eq.  (9.14)  may  be  rewritten  in  another 
form.  For  the  integration  in  h and  h’ , 


r r* 

<rj[ih'  cos  aecoj  Jh  » (9.17 ) 

-I  h' 

by  changing  the  order  of  Integration  and  then  inter- 
changing h and  h*,  the  result  is 


<rA'U' 


cas[Ka(h'-h)*ec6]<lh. 


'k 


(9.18) 


Referring  to  Fig.  (9.2),  expression  (9.17) 
represents  area  ~T~  expression  (9.18)  represents 
area  \/  . Since  cos  |K#(hf-  h)  sec  oj  is  symmetric 
or  even  about  (h*-  fa)  = 0,  the  total  area  may  be 


:»caa8S2B23K3E] 
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( • 


written  as 


I 


i M 

Ojfiit  a~n  eo»^0*-H)5*ceJ  dn. 


# ^ W M.  % 

\VmX9 l 


Then 


r 


jV^h-Ms-taJ  d h * Ojj’dJ/f  a^cosfc (h-h)sec6]<ih  (9.20) 

• f. 


°a“* 


-X  If 


where  the  limitation  on  (h*-  h)  has  been  removed. 

Thus  the  third  term  in  the  right  member  of 
Eq„  (9*14)  may  be  written  in  the  form 


-IfcpTTkf 


r fV  r*  . 

lajJMf  s*cOe  ctafalh^tee cejje  (9, 

J«JL  JD 


21) 


With  these  reductions*  the  force  R.^  becomes 


Rx  " X1  * X2  * X3  + x4  «*  X5*+  Xg'-i-  Xsn  (9.22) 


where  X^f  Xg,  X3#  X^*  are  given  by  the  equations 


O 


53 


following  Eq.  (9.12),  And  whore 


( 


x=»  - 

u 


- 16  Ttt>k* 


V * 16  f K, 


tr/ak'Jf'  crjliW!*«d36  • 1 ,ttC  °ce>* e7//d 

rt  n | *•  •* 

H-;,  •'.•t* 

Hr*-*1  r* 


Hr  Hf’-*1  f?  mfftf'bet'A 
$"  - - 32  TTf  K * <£&'#'  I We*  casfclMsec*]^ 


Por  the  force  R#,  mo  expressed  in  Eq.  (6.13), 
all  terms  are  real.  Hence  when  substituting  Hg  into 
Eq.  (6.13),  the  real  parts  are  to  be  taken.  The  final 
result  is 


R*  * Zl  ■*■  Z2  + Z3  +"  V * Z4  + Z5*  * Zsn  (9.23) 

where  Z Sg,  Z3  are  given  by  the  equations  rollow- 
ing  Eq.  (6.13)  and  where 


/'• 


Z4* 


wdhi 


♦ 
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z4" 


- 16  p 


K0J  jcr/jM'  6*6 
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• 52TPk; 
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r^iW 
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L-t  Jo 

There  is  no  corresponding  cancellation  of 
Z^1  and  Z4"  as  vas  ths  oase  for  Rx.  Also,  Z^1" 
cannot  be  reduced  to  the  integral  over  the  entire  dis- 
tribution A because  the  factor  sinj^fh*-  h)  sec  gJ 
is  not  symmetric  about  (h*-  h)  * 0 

By  following  a similar  procedure,  the  final  for 
form  of  the  moment  L is 

y 


Ly  * Mg  +M3  + «4  + lt5*  +■  M6»  + Mg" 

♦ *7  + *8  “9  + "li^iO^XO 


♦ «li+  *ll 
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where  H,2  M4  and  ^ given 

by  the  equations  follwolng  Sq*  (8*7)  and  where 


\u-6)<Tj!ilk'i4'  r«e£.36  e**^  6cjosfc(h'-4>)xcce}l6, 
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X,  Evaluation  of  Components  of  Rx 

In  the  foregoing  discussion,  the  expressions 
for  the  forces  Ry  and  Rz  , and  the  moment  Ly  act- 
ing on  the  toved  ship  have  been  formulated.  To  determine 
the  variation  of  these  forces  and  moments  on  the  towed 
ship,  the  various  Integrals  Involved  must  now  be 
evaluated.  The  evaluation  may  be  by  direct  numerical 
integration  for  a particular  case  or  the  Integrals  may 
be  evaluated  analytically  for  a general  case.  In 
either  case  the  amount  of  labor  involved  la  consider- 
able. In  this  and  the  following  sections,  the  analytical 
approach  will  be  used  so  that  later  numerical  calcula- 
tions will  be  applicable  to  any  particular  case  within 
the  limits  of  the  problem. 

The  present  section  deals  with  the  various 
components  of  the  force  Ry.  By  using  Eq.  (6.3),  the 
equations  for  the  source  densities  of  the  well-sided 
ships  A and  B,  the  first  component  of  Sx  becomes 


be2 

l* 


* mt 


l 

h'dh'df ' 


(10.1) 


-57 


Since  X^  is  an  odd  function  in  h*  with 
limits  of  integration  -l  to  l , X^  « 0.  This  is 
to  be  expected  because,  the  total  source  strength  is 
zero*  It  is  in  agreement  with  D* Alembert* s Pririelpl© 
which  states  that,  any  body  moving  uniformly  in  an 
unbounded  ideal  fluid  at  rest  at  infinity,  experiences 
no  resistance* 

The  second  component  of  Rx  is 


It  represents  the  resultant  horizontal  force  on  the 
distribution  A due  to  its  own  image  system  excluding 
that  part  of  the  image  system  which  trails  aft  to 
infinity*  To  show  that  this  component  is  also  equal 
to  zero,  consider  the  function 


h*h  (h»-  h)  # 

[(h»-  h)2f  (f  + f )2]^ 


f(h»,b) 


(10o3) 
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By  interchanging  h*  and  h 


f(h,h») 


hh»  (h  - h») 

[(h  - h»)2+  (f  + f»)2]** 


or 


f(h,h») 


h»h  (h»-  h) 

[(h»-  h)2-h  (f  «-  f»)2]% 


Thus 


f(h»,h)  « - f (h,h» ) 

Hence  f(h*,h)  le  skew  synmetric  about  \h*  - h)  = 
and  for  the  limits  -2  to  l for  both  h*  and  h, 

= 0 
A# 


The  third  and  fourth  components  of  Rx 
represent  the  resultant  horizontal  force  on  the  dia 

0 


•**«•■#***•  . w 


(10o4) 


(10o5) 

0 


(10o6) 
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trlbutlon  A due  to  the  distribution  B and  its  Image 
system,  excluding  the  trailing  distribution.  Because 
of  the  presence  of  Borne  common  terms  in  the  process  of 
integration,  it  is  convenient  to  ov&lunte  these  two 
components  together.  With  the  definitions  of  the  source 
densities  given  lay  Eq,  (6.3) 


X,  + X, 


4 <>!»*<* 
7 Tl* 


rd/J 


fc'dh'jn 


(10.7) 


’L-t 


where  «<,  and  d2  are  given  by  Eq,  (6.11).  Integrat- 
ing with  respect  to  f,  f*  and  h*  , the  resulting 
expression  after  some  reduction  is 


x3+  X4 


TT  l4 


rL+t 

(h-L)(l,+  ItUh 


L-i 


(10.8) 
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f (h - ?) [0>-i)V Ad1] ^ 


(h+i)+ni,ti)z+4j*]‘ 


j..'* 
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Th«  exact  integration  of  Eq.  (10.8)  with 
reapect  to  h is  a long  and  tedious  process,  the  final 
result  being  very  lengthy.  A process  of  numerical 
Integration  also  involves  a great  amount  of  labor  but 
is  more  suitable  beoause  the  result  is  more  concise. 
For  numerical  integration  Eq»  (10.8)  may  be  written  as 


■* 

>1  [-  8 (« UA,rAx+A,+A4) 

x 
z 

♦ * fc+'iKV  *1+  \o+  \i) 


Ifcf  b*cz 
TT 


\ 


(10.9) 


~ * ^rt)  J J 


where  ^ * (h  - L)/2l  and  o<  = L/2  l • The  quantity 
<A  is  the  ratio  of  the  distance  between  the  centers  of 
mass  of  the  two  ships  and  the  Bhlp  length.  With  the 
definition  fi  * d/2  l , the  ratio  of  the  draft  to  ship 
length,  the  terms  in  the  Integrand  are  defined  aa 


I r . -.2  ,7 

Ax=  + 4fi'\  , 


-Mxem ».  - 


. ..r 
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Z(S+  {fr  (<*•»$) 4(3*}  ‘ 


2 U+ 


A4«  ZfiUv 


2(3  {fcU+W + lj  %-  4(3*}  * I 


ZG***)-I-| 


As=  i:  [z^+VH]|[2^nH]*+4(3*JZ  , 


~~ir£z(f<+‘l)i‘,]j[2(dt>l)+tJ  +4&*J  , 
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A_-  2/3*1*, 


[z(<l+)l)-i ] + 4-4(3*}* 

[z(*+t)+tf  +{[t{d+})  + l]  +4<3*ji 


AS=  fM*)^  , 


r- {[z(^)+‘]*+ <Mzja  , 


A^q~  ~4 


4(5  + 


{[*(**■% )-l]  * HP*}  2 
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A1]l»  4<3jU 


4P+  [[2U+*l)l-lJZ+  IbP*}  * 


A12s  [[^*> -lj  + »*?*}  Z f 

A13  = “ if  (***)-m]  l60*j  , 


A14  = 


&+>?)-/]  + {[zUr*)-l]:+l6  0Z}2' 


[*(*+V+ 1]  t 


By  choosing  values  for  c*  and  /3  , Kqs 

(10.9)  may  now  be  evaluated.  The  procedure  is  to 

obtain  the  ordinates  of  the  Integrand  for  several 

values  of  ^ within  the  interval  - £ to  + $<,  Then 
» 

by  Simpson's  rule  the  approximate  value  of  the  integral 
may  be  found; 

This  force  is  analogous  to  the  mutual  at* 
traction  between  two  bodS-on_  As  is  to  be  exnected. 
preliminary  numerical  calculations  show  that  + x4 

is  9 very  small  quantity.  Individually,  however,  Xg 
and  X4  are  not  of  the  asms  order  of  magnitude  as 


srswMswaewiiKK 
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their  sum.  They  are  just  about  equal  to  each  other  in 
magnitude.  One  is  positive  while  the  other  la  nega- 
tive. Hence  their-  difference  is  very  small. 

The  component  Xg*  whioh  represents  the 
wave  resistance  of  ship  A as  if  alone  has  been  evalu- 
ated by  T.  H.  Havelock  for  a distribution  of  doublets 
instead  of  sources  and  sinks.  It  may  be  written  as 
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Vs-'iirP^ 


(P*+Q*)set*0d0 


(10.10) 


where 


cos 

Sin 


{h'-h)sec& 


.hi(WV\h 


(10.11) 


The  results  of  the  integration  of  this  com- 
ponent for  the  wall-sided  ship  described  in  section 
VI  are  summarised  in  the  Proc.  Roy.  Soc«*  Ser.  A, 

Vol.  108,  pp  582  - 591. 

The  component  Xg*  gives  the  resultant 
horizontal  force  on  the  distribution  A due  to  the 
forces  between  the  distribution  A and  the  trailing 
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system  of  sources  and  sinks  for  the  distribution  B. 
With  the  substitution  of  the  souroe  densities,  it  be- 
comes 


x * 
*6 


71* 


rd 


t fdfiL+l 


(h-L)dMf 
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s*t  ©<J6 


m(h-h)£Wd 
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where 
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. » K,3ecz6  Ism  >w({*4')l  - m coa  m ({+■£') 
jr(,M)s  _f k J 1 - (10.13) 

^ y ^ A 


•w  t*  *v$ec.  6 


and  (h*  - h)  < 0* 

The  integral 
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K6set‘d  «fj  - wcm  Mt  -Mxcoad 

c frxim  (10.14) 
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where  x > 0 has  been  evaluated  in  terms  of  the  two 
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variables  x and  f by  th®  Rational  Physical  Labora- 
tory (Teddlngton,  Middlesex,  England)  and  has  been 
tabulated  (10).  A numerical  method  of  evaluating  Bq* 

I 10*12)  is  possible  by  the  use  of  suoh  tabulated  re- 
sults* However*  it  will  not  be  pursued  here  as  the 
range  of  tabulated  values  is  not  sufficient  for  appli- 
cation in  Eq*  (10*12). 

Using  the  substitutions  v « m/  Ke  , 

’ |fa  (h  - h* ) » E,  K0(f  f 1 ) « D and  Integrating  with 

respect  to  ff  * Eq*  (10*12)  becomes 
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COS  V D -VH  cos  6 , 

e dv 

■v*+  cos4© 


(10ol7) 


The  limits  of  integration  Kdf  and  K,(f  * d) 
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are  to  be  substituted  for  the  variable  D.  By  using 
the  Laplace  Transform  of  cos  it  Eq.  (10.17)  may  now 
be  written  as 


M « — J fcos  VPcos  it  1 e HccS^dtf 

see1*)  *"  *-* 
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-gsecx6 
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By  writing 


c esi)X>  cosvg  s-  ^ v(g+p)  + cos 


and  taking  the  Laplace  Transform  with  respect  to 


r® 

^ F Hco**  + --Ji*£2JL---]e*SeC  6 &C 
2s#e  e|  [mWo^p)*  h We+te-o)2] 


The  derivative  of  Eq.  (10.19)  with  respect  to  D is 
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Then  by  substituting  Eq.  (10.19)  and  Eq.  (10.21)  Into 
eq.  (10.16) 
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dM 

dP 


H r I 

P)1  ~ 


(10-22) 


Using  the  transformation  A*  ? - D * Eq0  (10.22) 
becomes 


Msec  *6  - 4i)  = 
JP 


H 

sec© 


-P8«C*© 

e 


f®  A sec8© 

e 

iftoA+Aa 


-UD 


<lA . 


(1C023) 
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Then  from  Eq.  (10.16) 


-H 


r® 

e 


9 JL 
H1 

vLn 


H*!so*^0*Aa 


JA 


(10.24) 


Now  by  using  the  transformation  ten  0 » x end  Inte- 
grating by  parte 


i 


l 


l,  + la  “ 


tf<U 


0+k*)[h*+  (itxl)P*] 


(10.25) 


. 7 


-P(**X*) 

e <Jx 


v -A(i*x*> 

Jiif i\ 

Qt*+0+x*)Azj* 


By  letting  x be  e complex  variable  end 
choosing  a o on tour  consisting  of  the  upper  half  plane, 
the  first  term  In  Eq«,  (10.25)  may  be  evaluated  as 


p® 

HJx  __  vr 

r,.  » ,i 

<»*')[«*+ o+*z)vx] 

J© 

L ChN-0*)*J 

(10.26) 


It  can  be  shown  that  the  seoond  term  In 


*esw  aeis.u^aaagtsay  - r- 


C! 
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Eq»  (10x25)  is  of  the  order  l/H®  • For  practical 
purposes  the  magnitude  of  H is  auoh  that  this  term 
may  be  neglected.  Then  from  Bq.  (10. IS) 


«f  K.tV 

*14 


cl  iO.*t 


4 Vt 


r4_ 


H'dh'l  fa-LUh  [~-i — 

u 


ufuhtYk 

• • ^ f • • 


(10.27) 


M 


Now  by  completing  the  integration  with  respect  to 

h*  a t and  f*  and  using  ^«(h  - t*)/2l  , o<*L/2 1 

and  (3  = d /2i 


JL 

>!  *.»)+•*  V ®»4 

u 

+ f (^)(BU+Ba)-  8(^)(^,tBrt)]  df 


(10.28) 


whore 


Bg  * - Jz(a+*)tiJ  f |zft+*MiJ*+  40*}  * ? 


»— THU.. 


B3=  [z(*H?)-i]  , 


B5«  /S*jU 


£*(«h-4H]»  4pz]  * 


t {[*(*Hl)+l]V  4^]* 


Bg=  @*jbn. 


x •»  4 

Qf<A+if)->3  + { + H>&fZ 


sr_ 

“8 


s9=  {[2<4^>-!]% 

B10'  {[*t<^H]*+ «?*}*, 


wwwm»WMiK.( '.'  ilK88SEaag^-.^^^g  < * 
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Sn  « 


zU+\)-i 


®12  = ”£^h'I 


-j-  { [*W+*H 0%4^1} 1 


®13  = rl*> 


4&  + {fefttyi-Q  + «»£*}* 


2(*(+^)-| 


40  * (&<<+*)♦!]%- 

z(<A+fl‘M 


Then  by  following  the  same  procedure  as  for 
the  component  X^,  Xg*  may  be  evaluated* 

Finally*  the  component  Xgs  which  represents 
part  of  the  force  acting  on  ship  A due  to  the  wave 
Interference  caused  by  ship  B may  be  written  as 


Y n 
A6 


33PKBVc* 

TTV* 


h’jfc'jf' 

%B|> 

j,J.i 

oJl 

nLr» 


(h-LUhif 


Irl 


f?  - - r - 

sec.  6 e a&[*» (b~}»)  sec  6jj& 


(10.29) 
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-73- 


4- 

Integrating  with  respect  to  h,  f,  h*  and  f* 


X6  *" 


WPfcV" 


jjt(itfl9s*c*6&  cefyeoffasece) 


—o 


♦ (#*«»%  - *o/s)  COO  faLjtett)* 2 10  S»H  fcLtftc4) 


(10.30) 


+ (Kfl*u*6  - £A*%)ees(ig.x»cc6)  -2)0  a»*d*m  (l^L^eed)]  d$ 


where  L-,  « L ->  2/  and  Lg  * I.  t*  2 Z • In  complex 
notation  X6"  may  be  written  as 


- *'*«*  [,-i^AJlki;VU,aWe)ca^L.«ce) 


i7.  i i _£  r”  /wi 

r nmi  f«rb  «•  (J»»  w IV 


•Tij  « J 


Seel  port  of  fcos*6(>i0Zjet&ti)**‘*tl'*set6J  ^<40  a 


r (IG.ol; 


MNSSRmmMesK4eeaeMMMiMi 


'&s*mtei&ssu 
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4- 


To  evaluate  Xg"„  eonalder  the  fix*at  inte= 


gral, 


r-" 1 «/l§ * {£i.s*cs>  4$ . 


Ma  *a  \ 
\ XV»  VC  / 


Letting  p * 2 K6i  » Ke  d * & p where  * dj2l  and 
using  the  transformation 


aeo  8 « 1 + t/p 


(10*35) 


with 


d6 


it 


(,+ £)(>  + £)* 


(10.34) 


the  integral  may  be  written  aa 


# 


F - P..I  _f  fV^P 

'a~  »*tM*  pM  *.  v*  ■ ■■  .. 

1 


**  (l+aj)i 


X[(-V 


(10*35) 


T 


jRMMHfcS£suauiQCT;%-4*£4S^^  ^sss***. 
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where  t-4  * L/£  l . Each  term  In  the  expanded  Inte- 

grand Is  of  the  form 


^2*  * • 


.f* 

r 1 - , 4 

I •*("$>* 

- O'+f)4  p*0+f) 

i " 


(10.56) 


where  Y takes  on  the  values  3,  2 , ip  « 

A close  approximation  to  eq.  (10,36)  may  he 
obtained  by  expanding 


In  ascending  powers  of  1/p  and  integrating  terra  by 
term.  The  result  of  the  expansion  is 


For  practical  values  of  p « ZKtl  , the  expansion  to 


1 


pifiirrmrrr - ~iiinwWTiiirwMlMaiii>iFrBl  ~ 


1/p®  la  sufficient.  With  this  expansion,  each  term 
in  Eq.  (10.36)  Is  of  the  form 


4 


^r-i  -Or-U)t 


(10.38) 


where  r « 0,  1,  2,  3,  ....  • 

To  integrate  these  terms,  known  relation- 
ship (11) 


f 


-at  t“l  -jt 

« t <lt  * • p(z) 


(10.39) 


where  the  real  parts  of  x and  z are  positive  and 
■p  (z)  is  the  gamma  function,  is  used.  In  the  present 
notation  Eq.  (10.39)  may  be  written  as 


f 


Vi-drAW* 

t e 


o; 


r("i)*r*4«,t"MU  tio.40) 


where 


(*V •**)“*  , 

0*  6 re  ten  «y>  ► 


(10.41) 


HKF  £ * V£*-  ♦ 
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Uaing  Eq.  (10*40)  and  the  expansion  given  by 
expression  (10*37) , Eq*  (10*35)  now  beooraos 


JX*Feil  frt  J [$,  (*)  ^Q,**)]  (10.42) 


where 


<n<  ?) 


"s 


6 


fiOMjT  ».  li£W  v*r^.^ie 

[307*  * * !«:-i  7 4 * 


f (10.43) 


a.  97X7ryI4  |ld  isn8frf_fid 
10*4  r*  * f 4 e 


ss 

z 


% * »• 


/ 


with  ^*0,  2 ^ In  turn,  and  S end  6 

being  defined  by  Eq.  (10*41) 


Now  consider  the  second  integral  In  Eq. 


%* 


( 10*31)  and  let  It  be  denoted  by  Jg. 


J * R&ai  port  of 


jf[, -«*"**•]* 


(10.44) 


Using  the  same  substitutions  as  were  used  for  J-,, 


6 


_k  -jup 

Jo  - Pool  port  of  -LS 

4<2p>* 


X 


J.  t *(.♦«>*)* 


| (10.45) 


whore  </.' 9 L^/Z  l . Following  the'  same  procedure 

J2  * ?e«l  port  of  ^ ~y£ (10.48) 


where 
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A 


o 2(r)  - A~*lat4iA‘t,*]±. 

+[£**•'***  tgW.  (eg-*ir)si;f>‘ 


i£3£r*,t- 

SO JZ  * * e 


«*• 


»4otryvf-»t» 
1624  T ® * 


V (10.47) 


with  V = 0,  2/S  # 4^3  in  turn  and 

© “ arotan  »</  Z'. 

Similarly  for  tha  third  lntagral  in  Eq.  (10*31) 
denoted  by  Jg  , 


Jg » ?ea|  jm-t  efj  £|  - e * j 

X jc«sc$  (| secd  + t)***  K*Lz**c®J  J0 


\ 

’ (10.48) 


J 
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tha  result  Is 


?u|  ?art  ^ IT* [QaW-*c%3 fy)  (10.49) 


* 


»hsre  flfsLg /2l  , andQgty)  lathe  aaiae  as  Qg  ( y ) 

with  the  negative  sign  on  the  exponential  factor 
changed  to  a positive  sign.  The  definitions  for  S 
and  0 are  S*D'V(A"),J  ^and  0 • arc  tan  d</X  • 

Finally  by  substituting  J0„  and  J, 

• ' u 

into  Eq.  (10.31) 


X6  a~ 


32jj£c* 

ir^K* 


(10.50) 
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%’  XI « Evaluation  of  Components  of  R£ 

The  first  vertical  component  Z^#  the  resul- 
tant force  on  the  distribution  A due  to  the  forces  be- 
tween the  distribution  A and  its  image  system*  exclud- 
ing the  trailing  system  of  sources  and  sinks*  may  be 
written  as 


/•  1 


kMk'Jf' 


k0,  *’},«}{ 


(11.1) 


where 


f tf 


^ ~ [(h»-  h)2  *t  (f  + f»)2]^ 


Integrating  with  respect  to  ff»  f and  h* 


- 


4Pb8c* 
T ri* 


(11,2) 


I 
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* 


where 


c-(h)  • i f [(k-lAr] *- - [u-o‘r ■><)•]* 
t [miV-u*]*! + (tie)  ^otir*(*^) 


-cck-'(blij-i  ^"(g) +±C,tk-(lg)] 

■ aJh  f-tfink  t 


-Club' 


Ytf)} 


> (11.3) 


As  was  done  previously  for  the  component 
X3  t X4  , the  integration  of  Eq.  (11.2)  with  respect 
to  h may  conveniently  bo  done  by  numerical  integra- 
tion. By  letting  T « h/2  l with  /&  » d/2  t , 
Eq.  (11.2)  may  be  written  aa 


mmmwm 


,C  sri  - ■■ 


'.--Jr.- 


».  * - •vi  iWSr- 


-85- 


Z,  - - [p1+  p,+  p,  + v4~] 


+■  r0+-»r*)  [v,  t v+t  d7  t p8] 


V {11»4} 


8 (9 T*  {” J> . -*•  B + D b ll  Jt 

L *1  f ' *»  ■ -iij  r - • 


shsrs 


D3*~  [(*T-i)*-r  I4|»*]*  , 


D4*  [(*T+I)l+  U>0*]*  , 


Dg*  jt»v  ****  Hft'MV'fg*,]*  ^ 
(zr-i) 


n 


rsss 


*94* 


p.  - -u  jgJ*  , 


(VT-I ) 


8 


« j-tn.  4<g  + [jfrT»i)*+ 

z frr-H) 


De  * - 1*  ( (* r-i)  + [jir-> )%  4P2]  * j*  , 
Dio  = U ^*T+t ) + [pr+i  )S-  a?*]  * } , 

Dll®  /*,  |(zr-i)  + [(«•-!)**  *?*]*}  > 
D^g=  Z*v  |(zTH)  tjj[zT+l)V  l(>^]*J  * 


By  choosing  values  for  /£  and  using 
Simpson's  rule,  may  now  he  evaluated. 


T 


a5WJj5S*sr  r 
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%■  The  second  vertical  component,  Z2  , the 

resultant  vertical  force  on  the  distribution  A due  to 
the  forces  between  the  distribution  A and  the  Image 
system,  excluding  the  trailing  system,  for  the  distri- 
bution B , may  be  written  as 


wwr 

-i 


| (b< 
•K-i 


(11.5) 


where 


t 


f + f » 

[(h*-  h)2  + (f  + f»)2]* 


This  component  may  also  be  evaluated  In  the 

same  manner  as  the  horisontal  component  X_  + x.  • 

3 4 

Integrating  with  respect  to  f,  f and  h* 


Z 


2 


irM 


pi+i 

(H-l)  C(h)Jh 


Jlri 


(11.6) 


where  C(h)  is  given  by  Eq.  (11.5).  Then  again,  with 


r»T5®W*S 

1 
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f * (h  - 1*)/  tl  , 4x  h/2  l t (3*  &/21 


• I i»»  ? ® 

tor»*-c 


2 TT 


'Lx 


(11 .7) 


whore 


El*  » 


*2  * ~ {[2(4+*t)+l] \ , 


8, 


- fc(4+4H] 


% W**}  ^ 


B4*  , 


b5=  k 


zft  * 0X+'iP2}* 


B6*-V  !**♦ 

I a («***}  + 1 


T 


^■^MwisUwt f.-  «<****.*»,•..»  HBM&.  -.-if- 
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* 


B7=f  u \ A?+  + \l£± 

2{<A+H)-\ 


*>  4?  I 


z ^ 


K9*-iv 


E10=  ^ 


40  » f iT%  »*<9*)* 


[2^)t|]  t {[Z(4t^)tl]2t40i}t 


fat*)**]  + 


Again,  by  ohooalng  values  for  ol  and  /£  , 
and  by  using  Simpson1*  rule,  tbs  approximate  Talus  or 
the  integral  may  be  found* 

The  third  vertical  component,  the  resultant 
vertical  force  on  the  distribution  A due  to  the  forces 
between  the  sources  and  sinks  in  the  distribution  A 
and  B , is 


i z 


4?l\ 

Tl*» 


4 

4 * 

J0. 

•i 

(h-t) 


(11*8) 


38 


%P'  where 


* [(h*-  h)2*  (ft*  fM2]^ 

That  this  c ompoaont  is  equal  to  zero  le  not 
evident  from  Eq.  (11.8).  By  integrating  with  respeot 
to  f*  and  ft  Z3  may  be  shown  equal  to  zero.  The 
first  integration  with  respeot  to  f»  gives 


i 

where 


Z.  a 


TTZ41 


rX 


Wdfc’j 

J.J 


(h“t)  F iKdf 


(11.9) 


L-X 


p ^ ! ! (lioio) 

with  H*  ser  (h  • h* ) . 

Then  by  integrating  with  respeot  to  f , it 
may  be  shown  that 


'3 


0. 


(11.11) 
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turn* 
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%*'  The  component©  Z • t Z * . Z_r  may  do  evalu- 

ated  in  th«  same  manner  aa  Xg*«  By  following  the 
same  procedure.  It  can  be  shown  that  these  components 
are  of  the  order  i/*H2  • At  in  the  case  of  the  horizon- 
tal component  Xj»  , these  small  terms  will  be  neglec- 
ted* 

The  components  and  Zg*  oan  be 

evaluated  In  the  same  manner  as  * The  component 

Z^x  represents  one  part  of  the  resultant  vertical 
force  on  the  distribution  A due  to  the  forces  between 
the  distribution  A and  the  trailing  system  of  sources 
and  sinks  behind  the  distribution  A*  Substituting 
the  definitions  for  the  source  densities  it  may  be 

t 

written  as 


m 


3Z  PK.Vez 

TTi-f 


r) 

h'jK'ar 

•* 

• 

h (11*12) 


1 e dm 
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Integrating  with  respect  to  f,  f*#  and  h' 


# 


Z4m„ 


- ?a^yT  [ G-  ***^Jjj^-fl.*d$*"(F***) 


r (ii.i3) 


* co/d  C09(p*€C&)  + ■^•COidJd© 


where  p * 2*^  and  /3p  « . 

In  complex  notation.  Eg,  (11,13)  aav  be 

written  as 


Z 

4 


izdri^c8- 


l»wj.  pert  *f  j^*©*4©  +yr]  e **aec* 


r (11*14) 


+ 


CCtf  © 


J& 


J 


beoomee 


Using  the  transformation  (10,33) , Kq.  (11,14) 


(11.15) 


«.Srfc*  .,^-,^^|^-,..|-..j..-1  '-jfjiplmwi — r,T"  ’"-"‘TrT* — Tt* . 


where 


r® 

‘ill 

[<•«-  fl 

^0 

**0*f)‘ 

'b& 

l4 

(11.16) 


When  expanded,  eaoh  term  in  Eq.  (11.16)  la 


of  the  form 


-v 


+t)+¥T  » 


where  Y takes  on  the  value  O f 
How  by  expanding 

.-*!>*)  »£)‘ 

O'trO'i-X' 

in  ascending  powers  of  (1/p)  and  using  Eq.  (10.40). 
Eq„  (11.16)  may  be  written  as 


*3  * T*  r J - (*{3 ) V «*?Qa  bp)]  (11.17) 
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where 


*4< 


t > 


A 

1 e + 


ZS 


| 


I5f|-e 
tA 


l A'**. 


l fzi 


sgfOAt'*+^ 


*(iXr 


^ » * » 


with 

S *(**+»)"*  1 


0 - arc tan 


X 
* - 


How  returning  to  Bq.  (11.15) , 


£2uk*^i*s  •sdkmisui miwmam.  s?®^>*?ais«aMCSiS3^  •*■•*.•■* 


. ?4  ajfiSK&S 


S3 


* 


cos  Q JO 


(11.18) 


By  expanding  Eq^  (11.18) 


x4-|  u*6i6  - *| 

f^-alqd  0etli 

•I* 


co* 


(11.20) 


The  first  term  la  equal  to  one.  The  lest  two  terms 
ere  of  the  form 


eomBdO 


where  £ equal*  K0  d end  8 K*#  d. 

Bj  using  the  trensf  onset  ion  aee29  « It  t2 


V 


S « 


-I* 
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.2  _ r*2 

0«-t*)*e  it. 


(11.81) 
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>.V"  .1  wSWtfWWrfWt: 
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♦ 


Eq.  (11.21)  la  in  * foi*m  ouch  that  the  confluent  hyper- 
geometrio  function  (12) 


;^r " _***■ 

®k,m  f 5 * — (It*)  (itt*)  * * (x4)Ut  (11.22) 

r(* -•<♦**.) 


may  be  used.  When  k = m * 


W 


* 1 O-rffV** 


(11.25) 


i 


Henoe 


V 3<fJ  <«♦**)  * * R IV  (f  )1T*.  (11.24) 

*»  * 

Then 


I4  = hZW±  (K.i)TT*+  W,  , (11.25) 

*>  a *X»“t 

For  large  values  of  /£{  when  |«**g  < TT-^OT^ 
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the  asymptotic  expansion  of  a(  ? ) i®  given  by  the 
formula  (12) 


wk  _(  S 3) 

* C *»««  n}  g»  j 


J 

(11.26) 


When  k « a » -J  , Eq0  (11.26)  reduces  to 


-&*-*  * L d|  «$*  3M|*  J * (11.27) 


The  expansion  of  the  confluent  hypergeometric 
function  fop  all  values  of  £ such  that  | arg  f | <.  TT  ; 

and,  for  values  of  arg  £ ouch  that  JT  ^ J arg  £|  < -JiT 

is  (12) 


w, 

"k 


f < ) s — i 
’ XTTi 


>•» 


+i) 


(11.28) 


Whom  k e s s -* 


M.-k(  5 > 


*«>» 

r(«)r(-*+-t)r(-*+^ 


(11.29) 


-®  i 


T 


acr1* 


The  integral  may  be  evaluated  by  calculating 
the  residues  at  the  simple  pole  eh en  i And  at 

the  double  poles  when  a = 5/8*  5/2,  7/2,  . 

The  calculation  of  the  residue  at  the  simple  pole  is 
simple  enough  but  the  double  poles  pose  a more  diffi- 
cult problem  because  the  residues  Involve  logarlthmlo 

• 

terms*  By  using  a known  result,  Eq*  (11*20)  may  .be 
evaluated  without  calculating  the  required  residues* 
When  k = m = —5/2  Havelook  gives  the  ... 

result  (4) 


% i<  f > 

-f»-f 


jg.  iV1*'*' 


(ii 


where  jr  Is  known  ais  Euler's  or  Masoheronl's  constant 
and  la  equal  to  0*5772157  *..*  • ( As  a point  of 

interest,  the  value  of  J has  been  calculated  by 

J.  C,  Adams  to  260  plaoes  of  decimals*  ) (12) 

» 

Prom  Eq.  (11*50)  Havelock  obtains  the  result 


I 


5\ 

«3  de 


ill  r- 
3072 


r di.31) 


at? 

X048 


i*+ 


if 

1024 


A 


♦ 


By  differentiating  Eq,  (11.31)  twloe  with  reapect- to 
J , the  required  Integral.  Eq.  (11.20).  may  be  ob- 
tained without  calculating  the  residues  at  the  simple 
and  double  polea.  The  result  la 


ef 


cos  Be 


r (U.32) 


By  substituting  the  appropriate  expression 
into  Eq.  (11.19).  1^  may  now  be  evaluated. 

Another  method  of  evaluating  Eq.  (11.21)  is 
possible.  Integrals  of  this  form  can  be  evaluated  In 
terms  of  Bessel  functions  of  the  seoond  kind  of  orders 
aero  and  one.  Consider  the  integral  (13) 


« «•** 
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\ -***  *H±.' 


*j  «**( (itrt^t-  # ^(itt*}* 


(Xl«33) 


when*  n la  an  interger. 


This  say  be  reduced  to  tha  general  for® 


t * SL±£  L . .ga,“  A - 
2nrl  / 2e  -1  2|“  2n-  3 


Starting  with  n • 0,  It  oan  be  sean  that  if  and 

L _ ara  known,  all  tha  other  L’a  may  be  determined. 
•*1 

The  definitions  for  and  ara 


Ll" 


=39- 


To  evaluate  Eq,  (11,31),  from  Eq,  (11,34)  it  nay  be 
written  that 


I (11*37) 


By  applying  the  transformation  t ainh 
to  the  integrals  and  , they  can  be  trans- 

forma  d into  known  express ion a for  modified  Bessel 
Functions  of  the  second  kind  Kq  and  K*  (14) . 


\ (11.38) 


Since  the  functions  K0  and  are  tabulated  (14), 

L„3  may  now  be  found  from  Eq,  (11,37),  Henee 
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« 


-f 

e 


.1  _*#* 
% " Jt 


(11*39) 


The  letter  method  of  evaluating  1^*  appears 
more  oonolse  and  will  be  pursued  here*  When  $ = #rod 


V - f «W«  ^ (¥)] 


When  $ - * K*d 


V * M e 


[K(K.J)-Jr.<K.d)]  . 


Thus*  from  Eq.  (11.19) 


I4«  1 


(11.40) 


(11.41) 


(11.42) 
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I 


* 


How  by  substituting  I3  and  I4  , Eq.  (11.17) 
and  Eq.  (II*  \2)  roe; actively,  into  Eq.  (11.15)  the  com- 
ponent Zj"1  may  be  evaluated. 

The  vertical  component,  Zg"  t one  part  of 
the  resultant  vertical  force  on  the  distribution  A due 
to  the  forces  between  the  trailing  system  of  souroes 
and  sinks  behind  the  distribution  B and  the  dlstribu* 
tlon  A Itself,  may  be  written  as 

Integrating  with  respect  to  f* , f,  h’  and 
h and  writing  in  oomples  form 


(11.45) 


».  x 


Z5«a- 


32  Pb 


lZra>%+aD!t4&Jsm  fa L *ec6) 


+ tmog.  port  of  fcl *zt.Q  * l)*- (-£***"'***)  J 

+ litM^.port  of  jjcs*3  (k’J •*£*•* 


Je 


(11.44) 


a 


»a*wftxx»i wnw*  > —warn*  wawegg-aatM 


where 


*10?- 


t 


where  p » 2 K*l  , /Sp»*#d,  Lj  » L • 2 { and  11- 
How  by  using  the  transformation  given  by 
Eq.  (10.53)  and  using  Eq=  (10;40) 


» t.  *.  o 

9.  — ’ - 


2." « . sj±)lI  rSJ  +*+*1 

O l™6  * J6  * J7j 


(11.45) 


m 


$ 


with  the  following  definitions  for  J5»  J-  and  J-7- 


J5  -Iwj.f.rt  .f  (11.46) 


«s<  r ) - f 

♦ fef  8f*A*“.gjr 

"[«f  A*i*-«Trt,-<to*  «T ' «V“ 

* f vA**t  - 5 V'Jj,*  • • • • 


with  $ **  ( + o(  and  0 m arc  tan  ©</  Jf 


* 


Sg  a -<K*«rae  • ■ « rj; 


«■**:*.•  sior.T-  km>. 


C** 
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J6  - t mag.  part  fi)]  (11.47) 


4fcf) 


where 


v r » - }V11*.[^rA-»Vft»;^-4isV*w]a- 

+ [jtTAVK  Vfri  i;i  •’*,  (tg.  I 


[««**>$ 


m?)  v »)<*"** 


*»  -fid  1 | 

-t5tf  Jp*"" 


-1 

with  $*  jVV and  8 * arotan  d*/  Y • 


m 


^Imog.port  of  Si^j  T* [<),(») (11.43) 


■ >,  - -a 


where  Q^(  X ) la  the  same  as  Qg(  Y ) exeept  that  the 
sign  of  the  exponent  of  the  exponential  faetora  ia. 

__x 

positive  Instead  sf  negative)  and  where  f * [Y*+  z 

0 * arotan  / Y « 


.wmbu«*j(  jaua  Miriii  n I— f 1 HHti  U HU"  T m> i.h'i  r~~i  in  friTT~'~'  - f*JP. Aisrfltfu  i*.r 


,r  - * * 


■-  .vwcm.>WM»nK 
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XII.  Evaluation  of  Components  of  Ly 


Sines  the  integration  or  the  various  moment 
components  fellers  the  same  procedures  as  were  used 
for  the  force  components,  the  detailed  steps  will  not 
be  indicated  in  this  section*  Onl y the  final  results 
of  the  significant  components.  Kg*,  Mgn,  K10  and 
Mu",  will  be  listed*  The  other  components  may  be 
evaluated  by  the  procedures  indicated  in  the  previous 
sections*  The  components  Mg'  and  Mg"  represent 
moments  due  to  horizontal  foroes  while  and 

are  moments  due  to  vertloal  forces* 

From  Section  IX,  the  component  Mgt  may  be 

written  as 


V* 


Vc* 
Tl « 


t 


A 


(12*1) 


We  r ' x 

set©*  ce.*|Z#(h-l»)s«<©Jd0 


By  completing  the  integration  with  respoct  to  f,  f. 


# 


is&nmtt&iMMmifli?* WWK*""* 


.»lv; 


h,  h'  and  0,  Mft*  may  be  written  ao  the  sum  of  throe 
parts. 


V * • 

?hs  txpriatloni  for  *«'•  “ss’*  and  are 

"si’  * - {£* * 

®*  • 

I t,4|  1 r *l^»^  fas  /a\_  .r^(n 

l'e*‘  F*1^  • • * ^'vf  **  *y~rqj  5 

sr  » _ l4Ptffc*#  f P*.  (4 t.-PP.  #^  \ , r^P. 

"52  {t  * £ re  L-y  w 

r^ir^^a^p)  + *e  Wt7^a^p) 

♦SW  f»rt  •»  Jg^T1  ftfrj-w't.pf  ] , 

-53*  - w 

- e*,PLs.(*?f>) -*•  ' **FLy(z?r ) 


(12.2) 


(12,3) 


t ?«!  fart <f  [e^ - **%<**)] | > 
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m 


t 


where 


L-b  (?)  * £ . 

L-7  (?)  •%'*[!(-!*■$*.(£) 

^*■•9  -^+3^.(4) 

+ *(-r ^(|)j  , 


and 


<y  y >-  sKil°]  + 

-&iS***i6-4  s V~  ] -L  + L iZZLf  v*r  C * i0 


1—2+1-  i^22 

J Fz  [ 307 z 


- sar  r*s  *,  % c §f  j j-*,  zaf  r ) s t 1 « 

rCZ4  * " ‘ * 


\ V4  jo*-?  / 


(^,V-iBU-€)  SV'\  /js  rt  a«i  \ ;ff'« 
\ /«>  50 7i  / l * * V c 


+ f S 


$e*i€>  7 JL 

J P3 


• • — 


J 


(12.4) 


(12o5) 


In  Ctg(  y),  i*  { IfN*  1)*^  and  0 •»  arc  tan  »/y  . 
Hj(  Y)  la  obtained  from  E^*  (10.47)  as  explained  after 
B<j*  (10.49).  However,  when  Q^(  Y)  la  used  In  M52* 
and  Mw*  . 5 = { rH  1>"*  and  3 * arc  tan  1 /Y  . The 

£'s  in  the  L(  J ) function*  take  on  the  values  p 
and  as  Indicated  in  Sq.  (12.3).  The  functions 
K0<  3f  /2)  and  K^(  $ /2)  are  modified  Beseel  functions 
of  the  second  kind. 

The  component  Mg"  may  be  written  as 


V 


»3PK.*bV- 

Tl< 


r^r1 


••w 


oLtl 

(K-t)JWf 


L.-1 


5«  0e 


[kjt'-hkKo]  de 


}■  (12.6) 


J 


or 


62 


M " 
63 


(12.7) 


where 


w 


>Mwa»iaf«wen 
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-rW  part  af  j^E**  JcjSi(0)-;e^ri?|0(^)+ 


+fe.l  rrt  .f  -«A}{  ftf> ) (ip)]  | 


I » * S£ 
62  Tfc 


i-Vfe  r OT 

w L 8,i 


* J2  * JS  I » 


■] 


u 32PI* 

63  TCK 


^pr  {«•>•  m [;»Qlw-;w<;,<ie)] 


yort  of  jfo 

•+  ^>}  f»rt  »T  | • 


For  the  Q( ¥ ) functions  in  Eq«  (13*8)  — 

Y)  is  given  by  Eq.  (10.42);  / ) is  given  by 

Eq.  (10.47);  end  Y ) is  given  by  Eq.  (10.47) 


•i~rr-inr  ~ - iiiunwnwupm  '/  .r  i si  n — nf  n~n  irnr  rrm — 


(12.8) 


-no- 


except  for  a positive  instead  of  a negative  sign  on 
the  exponent  of  the  exponential  factors*  and 
$ - (<*")*]  0»  arctan  <*? / Y • The  definition  for 

%(Y)  is 


v r >-  [f  «V*.  f *«.«••]  - 


4S’e 


i *i* 


M 


"i£y3S  *e 


^ i© 


+ 23fifyaf^^  222!£Vff*.*id 

!S2e  - * - • !fl24  ' ' - 


iCZ-4 


(12.9) 


+ KH4^  sM‘%.  |rsf.fi ‘V§ 


947Z 


where  $*(**+  <**  5 ^ and  0 « arc  tan  o</  If . Q^0(  y) 
le  the  same  as  Qg(  given  by  £q.  (12.5)*  except  for 
a negative  instead  of  a positive  exponent  on  the  expo- 
nential factors.  Also,  the  definitions  for  S and  6 
In  Qg(  Y ) are  S = i/Vd*/*]  * and  6 * arctan  <*/  y . Pur 
Ql0(  y ) the  definitions  are  f*  0fS-(4')xJ  **  and 
0 * arctan  <A.'/Y  . 


;^;«ap!»*W*SUWl 


For  the  three  J functions;  J^  ie  given 
by  Eq.  (10.42),  Jo  Is  given  by  Eq.  (10.46),  and  J, 
is  given  by  Eq.  (10.49). 

The  component  K^1  may  be  written  as 


»» 


>2  pCfcV 


10  IT  l* 


»4 

rZ 

X 

J 

** 

(itfW 

hdhJt 

i J 

o 

h • 

-*,(f*'W*  r , n 
0 « iitijjc/h-KjjecejclO 


(12.10) 


and  reduced  to 


*10  = •*  [<?»(*)- 

i-Vtolf art  of  [<?u(d)  - * e *^^(4(3)] 

+ t(T+f)-  "T « * (f  p) 

+2.  ez*?Lg.(2pp)  + *e~**Pt7  (app)| 


U<s«.u/ 


where 


* 
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♦ 


«u<  n-*  [£  «**•«■  § * V1*] 


tfiSf  &***"■>  -!3£|M'*+  iC  jM'1* 

J_*<*8  * |0M  ao48 


Jp-  l(,3*4 


t 2*s?si  H «*»•«*.**- 1 rjf«?is 

*384  -ttfcft  8 


4 0 9 0 


.u,M's]  jj 


0l3  nr ) = ft*.*'*-  [ft  lr  A***  g »*«***] -jr 


+ ,r»]£Sfe^6+ 

+ L«M1  **?*  +|0H*» 


^•»J-.ri22Sf  3 « ^i6 

J p*  |_«m. 


4§nrvV* 
sS?*5  * 


*l» 


«U£T  «*«**• 

5»12 


♦ 


IJI’TCj'.  i 4 «' 

« — 3 * 
on* 


* 


•»***i.rww«cs^r: 
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with  $■(  1)*^  and  © n&rotan  I /Y  % 

The  L function*  In  Eq*  (12*11)  are  given  by 
Eq.  (12*4)* 

Finally 


f 


Xp*« 


* 4 -*(f  ♦*')#«**• 


#i*i  r^(K'-^)*e<  ftj  JB 


(12*13) 


becomes 


V-  [J1S^14*  •'IS-  'l6—  *„] 


(12.14) 


where 


J13  -Ke*l  part  #f  T*  [^(h)-ae-^M  , ] 


J14*‘H^  of  (*&]  y 


(12*15) 


TT~ 


-114- 


t 


JX6  = ^ f- rt  V*  fyo)-Z«0FG/J*r) +**'%t(4fi)J  , 


rl7«K!Mfairt«f  f ‘ j^(o)-ze*f<^(**)  ♦ e *%„**)]  . J 


\ (12.15) 


The  Q{  Y ) functions  in  Eq„  (12.15)  are 
defined  as  follows: 


♦«•**“]* -[Jgp* 


1 fi* 


+ 2«f  rV*-4  **♦  *2# 

*446 


2^sV'\  *s *«*“♦  f Atu]^+....  I 


(12.16) 


where 


$«()T2t«<l)',“^  and  0 * arc  tan  oC/  X • 
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i4<  '<  > - <b*'**- 

'* 1 ' {'*]  f,  -[ *ar At ¥'*+  asHrtV*'*  . (12.17) 


:£"?•*  JL2c^"^,el  1 


+ «2S*VstH  isssfjl.- 

«m  *45-7*. 


where  4 and  8 ■ arc  tan  *><7  ^ » 

<*15  ( 2f)  Is  obtained  from  Q14(  if)  by  chang- 
ing che  sign  of  the  exponent  of  the  exponential,  factors 
and  by  using  [Y%  and  0- arctan«*V  JT  for 


$ and  & • 


<W  * ) - ■is#*- [irsf«fi%  gs V‘*J  J- 

♦ • (12.10) 


Jf1  L««sfi 


H 
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m 


4.  — — — a X # 4 

* *>H  ' ' C 


4*‘ !•» 


* 2£2£fi*e#i*-  J.ys'f.f16-  4<rsIJie  lx* 
1X288  4 4 J P* 


'Z 

where  S =^2+(^.,)iJ  and  6 *■  arc  tan  «</  f » 

«17<  r ) Is  the  same  as  Y)  except  for 

the  definitions  of  S and  6 in  which  sC/#  is  subati- 

. x 

tuted  for  <*'  , i.e.,  S'OrV^")1]  «d  9 = arctan  «*7  &f . 


r 


XIII.  Numerical  Calculations 

As  mentioned  in  the  introduction,  as  far  as 
the  effect  of  waves  is  concerned,  the  force  and  njOiSdnt 
experienced  by  the  trailing  ship  nay  be  niscaed  into 

Vt>-  £k  Ot  na«f  a • fVa  A «»•»»  #4  AL.  - * — — — * 

c 9 ~ v vuv  ouay  -a 

alone,  the  force  and  moment  due  to  mutual  interaction 
between  the  two  ships,  and  the  force  and  moment  due  to 
wave  int<  rf erence  • 

For  the  ship  when  alone,  the  horizontal  force, 
or  wave  resistance,  has  been  calculated  and  plotted  by 
Havelock  (4).  It  is  an  increasing  oscillating  function 
of  the  velocity.  As  the  draft-length  ratio  is  increased 
this  force  increases.  To  the  authors’  knowledge,  the 
vertical  force  and  moment,  due  to  waves,  on  the  ship 
when  alone  has  neither  been  calculated  nor  plotted. 

The  form  of  the  functions  involved  in  the  vertical  force 
and  moment  are  very  similar  to  those  m tne  expression 
for  the  horizontal  force.  Thus  it  is  to  be  expected 
that  the  vertical  force  and  moment  are  alec  Increasing 
oscillating  functions  of  the  velocity. 

The  force  and  moment  components  which  are  due 
to  mutual  interaction  vary  inversely  as  the  square  of 
the  distance  between  the  two  ships,  and  hence,  they  may 
be  expected  to  diminish  very  rapidly  with  the  distance 
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between  the  two  ships.  Pop  any  practical  distance 
between  the  two  ships*  these  components  can  be  neglected. 

Due  to  the  motivation  of  this  problem*  viz.* 
the  problem  of  the  effect  of  & leading  ship  on  a trail- 
ing ship;  of  major  interest  hors  arc  the  force  and  moment 
components  due  to  wave  Interference.  The  significant 
ones  are  designated  by  Xg%  2g"  and  Mg^",  Hgg" * 

Mgg" , W-n".  The  interpretation  of  these  various  com- 

ponents In  terms  of  the  distribution  of  sources  and  sinks 
which  replace  the  ship  Is  given  In  Section  VII  and  VIII. 
The  numerical  calculations  and  comments  are  confined  to 
these  force  and  moment  component*  = 

f*  The  parameters  Involved  In  the  numericel  cal- 
culations are:  (1)  the  ratio  of  the  distance  between 

ship  centers  and  ship-length*  oC  * L/2  l * and  the  corres- 
ponding  quantities  oC=  Lj/2  l * o<  - 1 and  <*  =L2/2  l = 

■sC.  f 2 i (2)  draft  ship-length  ratio  p=t  d/2*  » 

(3)  ratio  of  the  distance  from  the  origin  of  the  rec- 
tangular axes  to  the  center  of  mass  of  the  ship  and  the 
ship-length*  6 /2l  ; and  (4)  the  Inverse  of  the  square 
of  Proude*s  number,  p *•  l/F where  P « c/l/SgT' • The 
values  chosen  for  these  various  parameters  are  as 
follows: 


wr*  A. 
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f 

o(  * 2 through  7 
& « 0.1 
e /2Z  « 1/40 

p * 20 

Omitting  the  details.  Figures  (13*1)  through 
(13.7)  show  the  variation  of  the  various  components 
listed  above  with  respect  to  the  spacing  between  ships. 
The  curves  are  all  of  an  oscillatory  nature  and  are 
camped  very  slowly  with  respect  to  spacing  between  ships. 
The  frequency  per  ship-length  of  the  oscillations  in  the 
curves  remains  practically  constant  with  increasing  dis- 
f tance  between  ships.  It  is  surprising  that  even  for  a 

relatively  large  distance  of  six  ship-lengths,  where 
one  might  expect  negligible  wave  interference  effect, 
there  is  an  appreciable  value  for  the  separate  compo- 
nents. 

To  obtain  the  relative  magnitude  of  the  wave 
interference  off  sot-  on  th*?  trailing  ship,  consider  Xg1*  - 
From  Fig.  (13.1).  the  value  for  the  wave  interference 
effect  at  a ship  spacing  of  one  ship-length  Is  about 
0.00375 p b2g  l • At  a spacing  of  six  ship-lengths,  the 
value  decreases  to  about-  0.002  p b^g  l . The  wave-resis- 
tance. for  the  same  parameters,  for  the  trailing  ship 

& 
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as  if  alone  :i»  approximately  equal  to  0.005  /»b2g  l (4). 
Thus  this  component  of  the  wave  interference  effect  at 
one  and  six  ship-lengths  is  about  75 % and  40£,  respective- 
ly* of  the  wave-resistance  of  the  ship  as  if  alone*  Then, 
since  the  wave  Interference  effect  is  oscillatory,  by  a 
judicious  spacing  of  the  two  ships,  there  is  an  appreciable 
reduction  in  the  wave-resistance  of  the  trailing  ship* 

For  the  other  force  and  moment  components,  one  may  expeot 
similar  relative  magnitudes* 

At  the  speed  chosen  in  the  calculations,  all 
the  curves  oscillate  at  the  rate  of  about  three  cycles 

por  ohip—  leugtu*  Tlio  o^uiitiuuo  tuo  toTIOuo  forces 

and  moments  indicate  that  as  the  speed  Increases,  the 
rate  of  oscillation  of  the  curves  decreases.  Thus,  at 
a low  speed  the  wave  interference  effeot  on  the  trailing 
ship  may  ohange  three  times,  per  ship  length  spacing, 
from  a positive  to  a negative  value*  At  a higher  speed 
the  change  from  a positive  to  negative  value  may  occur 
only  once  per  ship-length  spacing  between  ships*  This 
would  indicate  thst  if  there  is  to  be  any  sustained 
saving  in  power  which  would  result  from  reduoed  resist 
tance,  it  would  more  likely  be  at  the  higher  speeds* 
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